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Abstract 

Montroll had considered a Verhulst evolution approach for introduc- 
ing a notion he called "social force" , to describe a jump in some economic 
output when a new technology or product outcompetes a previous one. In 
fact, MontrolFs adaptation of Verhulst equation is more like an economic 
field description than a "social force". The empirical Verhulst logistic 
function and the Gompertz double exponential law are used here in or- 
der to present an alternative view, within a similar mechanistic physics 
framework. As an example, a "social force" modifying the rate in the 
number of temples constructed by a religious movement, the Antoinist 
community, between 1910 and 1940 in Belgium is found and quantified. 
Practically, two temple inauguration regimes are seen to exist over differ- 
ent time spans, separated by a gap attributed to a specific "constraint", 
a taxation system, but allowing for a different, smooth, evolution rather 
than a jump. The impulse force duration is also emphasized as being 
better taken into account within the Gompertz framework. Moreover, 
a "social force" can be as here, attributed to a change in the limited 
need/capacity of some population, coupled to some external field, in ei- 
ther Verhulst or Gompertz equation, rather than resulting from already 
existing but competing goods as imagined by Montroll. 

PACS : 89.65.S, 89.75.Fb, 05.45. Tp, 89.75.k 

1 Introduction 

Malthus [2] had assumed that the growth of a biological population was strictly- 
proportional to the number of members in the population, dN/dt — r N, thereby 
obtaining an unreasonably excessive exponential growth N(t) ~ e rt . To correct 
for such a behavior, Verhulst attempted to find an evolution equation with a 
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limiting growth term [3J 0] . He empirically introduced a quadratic term, in the 
right hand side of Malthus equation, in order to mimic the food availability 
and/or the so called carrying capacity M of the country, 



He obtained the famous, sigmoid, logistic function. In contrast, Gompertz [5] 
formulated the decrease in the number of members of a population in terms of 
a first order differential {growth) equation of Malthusian form, but with an ex- 
ponentially decaying birth rate, r = r§e~ Kt . This results into another analytical 
sigmoid curve, though asymmetric, in contrast to the symmetric logistic map. 

In 1974, Montroll and Badger [5] had introduced, into social phenomena, 
quantitative approaches, as physicists should do along mechanics lines. Nowa- 
days, one would say that (i) "degrees of freedom" may interact with "external 
fields" , like a spin interacts with a magnetic field or a charge interacts with an 
electric field, and (ii) " degrees of freedom" may interact with each other within 
some cluster, e.g. like for a spin-spin interaction. 

It was observed, e.g., in [T], that deviations from the logistic law [21 H] were 
often associated with unusual intermittent events: wars, strikes, economic pan- 
ics, etc. Mathematically, the effects can be abstracted through an instantaneous 
impulse function, called a "social force", by Montroll Q], inserted into Verhulst 
equation. In so doing, some remarkable shift occurs in the time evolution of 
events, after some delay, - the shift size depending on the force/impulse strength. 
Montroll presented semi-log graphs, with such apparently drastic shifts, on the 
replacement of rail service by air service in intercity passenger travel, or the re- 
placement of sailing ships by steamships, - both rather "economic conditions", 
or an accelerating decline in agricultural work force, - an effect apparently due 
to telegrams from the White House. 

However, the notion/definition of some "social force", as in Montroll pQ, is 
rather loose. Nowadays, the notion of "social field" seems to be a more realistic 
concept [3 |8]. In fact, Montroll discussed industrial evolutionary processes, 
considered as a sequence of replacement or substitution processes. It is like a 
competition between two equilibrium states in a free energy-like description. 
In such a spirit, shifting the "thermodynamic state" toward one or another 
minimum, is rather due to the rise of an "economic (potential) barrier"; thus a 
wording like "economic force" could have rather been used. 

Incidentally, the analytical form of the force impulse as introduced by Mon- 
troll into Verhulst equation, see Sect. 3, appears rather ad hoc. Moreover, it 
only allows for a constant and identical growth rate on both sides of some time 
interval, - as that was the case in the data presented by Montroll. 

In view of such considerations, I was interested in finding whether one could 
get some "social data" on the evolution of some community, made of agents 
having a well defined " degree of freedom" , asserted to some external field, and 
whether one could observe some (keeping the same wording as Montroll) " social 
force" through appropriate data analysis. Moreover, to find smoother evolu- 
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tions and various growth rates should be a plus. This would request more high 
frequency data than yearly ones, as in [1]. 

An indirect measure of social behavior, like the construction and inaugu- 
ration/consecration of temples by a religious movement, seems reliable, more 
than, e.g., adhesion numbers as in [H3 OH [HI EH U: - relying on statis- 
tics based on interviews and surveys. Moreover the former data involves more 
immediate mixing of economy and social conditions than the number of adepts. 

As a practical example, the evolution of the number of temples, built in 
Belgium, by a community like the Antoinists [15] is reported in Sect. [2] and 
first studied through the most commonly accepted kinetic growth law, i.e. the 



Verhulst logistic function [3J; see Sect. 2.1 



Waiting for a change in taxation laws in Belgium, the construction of tem- 
ples stopped during a while, after the end of the first World War. When the 
constructions resumed, the growth rate apparently changed. Montroll's "for- 
malism" is thus inappropriate, on two counts at least. In fact, it will be seen, 
in Sect. 3, that one can introduce another analytic form of a "social force" into 
an evolution equation in order to describe a rate change, - with or without a 
gap in the evolution. 

Moreover, the Verhulst mapping is sometimes criticized as unrealistic. The 
Gompertz (human mortality) decay law [5] is often suggested to be the al- 
ternative, - tuning the parameters into their size growth value, as often done 
nowadays. Some quantitative analysis is thereby found in Sect. |2.2| along such 
a line of approach. Thus, quantitative measures, as expected within Montroll 
and Badger pioneering framework [6], are obtained. They are commented upon 
in SectlU 

In fact, Verhulst and Gompertz basic evolution equations are particular ex- 
amples of a more general first order differential equation [IB] which is discussed 
in the conclusions, Sect|4j with the idea of presenting the general addition of an 
external field-like term in such a framework. 

In brief, this study has a marked difference from that of Montroll. It is not 
a competition effect between two processes, that one is looking at here below, 
but rather an evolution due to an external field, - a taxation system here. Here, 
the "social force" is attributed to a change in the limited need/capacity, in 
Verhulst or Gompertz equation, rather than resulting from already acquired 
goods competing with new ones. Thus, it can be pretended that one can be 
presenting two different views of a " social force" , - as resulting from two different 
origins of the minima in some sort of thermodynamic free energy potential well 
description: a competitive interaction or the effect of an external field, leading 
to the choice of the (kinetic but equilibrium-like) state. 

Moreover, in the conclusion section, the mathematical form of both types of 
forces receives some sociological interpretation. This can be based on the change 
in the limited need/capacity of some population, coupled to (or influenced by) 
some external field, rather than resulting from already acquired but competing 
goods as in [1]. 
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Figure 1: Logistic fits, at low and high t, of the number of temples of the 
Antoinist Cult in Belgium as a function of the number of months since the 
consecration of the first temple on Aug. 15, 1913 in Jemcppe-sur-Meuse. The 
"size" of a year is of course 12 months. Each year is "centered" on July 01. 
N.B. "1910" on the figure is a misprint: it should read 1913. 




2 Data Set 



The Antoinist community exists for about more than a century |17[ I18j . has 
markedly grown, has expanded in various countries, but and is now apparently 
decaying. Although it might be of interest to consider the data for the whole 
"sect" on a world wide basis, only the 27 temples constructed in Belgium during 
the original growth phase of the cult are here considered. The date of every 
temple consecration has been extracted from the archives of the Antoinist cult 
and from [T7] and [TH] compilation and discussions. Most of the times, the exact 
day of the consecration is known, - twice, only the month is known. In such a 
case, I took, the 15-th day of the month as the consecration day. To be more 
precise about the exact day of the event, for both cases, has requested much 
time consuming research for this information through news media archives, - but 
without any success. Sometimes, some disagreement, between dates, was found. 
When in doubt, the dates in an Appendix of a 1934 book by Debouxhtay [17] 
were preferred, since I consider them as the most reliable ones, - due to some 
"proximity effect" of the book's author . 

In order to keep the same units, after calculating the number of days between 
two events, the number (m) of months (months) has been next rounded up to 
the nearest integer. 

Thereafter, the evolution in the number of temples was analyzed in terms of 
the Verhulst (logistic) law and the Gompertz (double exponential) growth law. 
In both cases, it readily appeared that two regimes had be investigated: (I) one 
between 1910 and 1919, for 16 temples; (II) another between 1923 and 1935, 
for 11 temples. For further reference note that the X15 value with 15 degrees of 
freedom and the Xio with 10 degrees of freedom have a critical value equal to 
24.996 and 18.307 respectively, for a 0.05-level test [TS]. 



2.1 Evolution study. Verhulst equation 

Let us take the (Verhulst) logistic function, with N/M first approxi- 

mation of a growing entity, i.e. the so called logistic map, a sigmoid curve, 

e r(t-t ) z 
Z ^ = Z °° 1 + e r(t-to) = 1 + e -r(t-to) ' ^ 

where z x is the upper limit of z as time tends to infinity, to is the position of 
the inflexion point and r is the supposedly constant growth rate. This way of 
expressing the logistic curve has the advantage that the initial measure, here 
zq = 1 at t = 0, is a rapidly fixed value for one of the three model parameters. 

The upper value z^ is a priori imposed to be an integer. The low t logistic 
fit of the number (N) of temples as a function of the number (to) of months 
(cumulated since the rise of the first temple, on Aug. 15, 1910) corresponds to 

N(m) = 24/(1 + e -0.03396*(m-80)) (3) 

while the fit in the upper regime corresponds to 

N(m) = 29/(1 + e -o.oi95*(m-i40)^ (4) 
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Figure 2: Gompertz double exponential law fit of the number of consecrated 
Antoinist temples in Belgium as a function of the number (m) of months (cu- 
mulated), in low and high t regimes since the consecration of the first temple 
on Aug. 15, 1913. The "size" of a year is of course 12 months. Each year is 
"centered" on July 01. N.B. "1910" on the figure is a misprint: it should read 
1913. 




time (months) 



Both data and fits, in the appropriate regimes, are combined and shown in Fig. 
1. For statistical completeness, let the x 2 values for the fits be reported: they 
are respectively 1.208 and 0.392, much below the critical 0.05 test value, recalled 
here above. 

It seems worth pointing out that the initial growth rate is about 0.03395, 
i.e. largely more than 3 temples every ten years, but is reduced to 0.0195 in the 
latest years, i.e. about 2 temples per year. Note that a unique logistic curve 
would give a value of the growth rate ~ 0.02355. The curve is not displayed 
since such a value does not fulfill the Jarque-Bera test. 

2.2 Evolution study. Gompertz equation 

Gompertz's law is also a 3 parameter expression, i.e. 

y(t)/y M = e~ be ^- rt \ (5) 

where Dmi b and r are positive constants. This corresponds to an exponentially 
decaying birth rate r in Malthus equation, i.e. r — roe~ Kt , pending ro and K 
being positive constants. 

The derived differential equation is commonly referred to as the Gompertz 
equation, i.e., 

(6) 

where k has mutatis mutandis the same meaning as M, the carrying capacity, 
in the Verhulst approach. 

Fits to Gompertz double exponential law can be searched for with different 
techniques H3 123 HI EH E2H |2B] , still imposing the amplitude y M to be an 
integer. It has occurred after much simulations that two distinct regimes must 
be considered, exactly as in the analysis along Verhulst approach: one at low 
time, i.e. during the initial growth of the Antoinist communities, and another 
at later (high) time, with a 4 year gap, between 1919 and 1923: 

N(m) = 23 e -.-«— >'«" (7 ) 

N(m) = 31 e -e-c~-"«»/™ (8) 

respectively, as shown on Fig. 2. Note that the upper (absolute) values of the 
possible number of temples which could be asymptotically expected, slightly 
differ in either Verhulst or Gompertz approach. 

Recall that the logistic law leads to expect ~ 24 temples at saturation, for 
the initial regime, but 29 temples at most, for the second regime. These values 
are 23 and 31 respectively for the Gompertz law fits. 

For statistical completeness, let the % 2 values for the fits be reported: they 
are respectively 1.480 and 0.249, much below the critical 0.05 test value again. 



dy 
dt 



r y log 
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Figure 3: Logistic variation (X/(l — X)) of the number of temples (X) in 
Belgium as a function of the number of months (cumulated from the raise of 
the first temple), in order to indicate the effect of a "social force" at time r, 
accelerating the process over a time span At a = (t a — t*) allowing to measure 
the strength a of the force impulse in the sense of Montroll. 
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3 Quantitative measure of "social forces" 



At first sight, the growth regimes do not seem to overlap much, to say the least. 
They occur on different sides of a time gap. Moreover, the rates of growth seem 
somewhat different in the successive regimes, indicating sequential rather than 
overlapping (or/and competitive) processes, as in Montroll report. 

It is worth to recall that Montroll argues that social evolutionary processes 
occur as a sequence of new ideas for old ones, inducing deviations from the 
classical logistic map associated with intermittent events. Montroll has argued 
that the most simple generalisation of Verhulst equation, in such a respect, 
goes when introducing a force impulse, F(t) — a d(t — r), in the r.h.s. of the 
kinetic Verhulst equation appropriately rewritten to emphasize the growth of 
the logarithm of X = zj z x , 

^"(1-4^^(1-4 (9) 

so that the dynamical equation for some instantaneously forced, at time r, 
evolution process X reads [27] 

dlr ^ X l=k{l-X)+a5{t- T ). (10) 

In other words, 

^ =kX(l-X) + aX 5(t-T). (11) 

In so doing, in the time regime after the withdrawal of the intermittent force, 
the evolutionary curves are parallel lines, on a semi-log plot: the unaccelerated 
one being above or below the latter depending whether the process is accelerated 
or deterred at time r. The impulse parameter a is obtained from the step-like 
shift [J to be 

a = k [1 - X(t)} At A (12) 

where At a = {Ia — if) measures the time which has been gained (if AtA > 0) 
in reaching some (arbitrary) A level, i.e. Xa/{1 — Xa)- 

Two such impulse effects can be seen in Fig. 3, which is an appropri- 
ate replot of Fig.l. An acceleration at t\ ~ 36 (months) and a deceleration 
after t 2 ~ 106 (months). One easily calculates that X(ji) = 0.14, since 
A(tx)/(1 - X{n)) = 0.16 and X(t 2 ) = 0.55, since A(t 2 )/(1 - X(t 2 ) = 1.23. 
To estimate fc, note that 125 months were required for X/(l — X) to be mul- 
tiplied by a factor 10 in evolving from 0.1 to 1.0 or from 1.0 to 10.0, so that 
k = (1/125) Zn(10) = 0.0184/month. Moreover, AtA = ~ At B = 40 months. 
Hence, the impulse strengths are: u\ = 0.103 and a 2 — —0.405. These are very 
reasonable orders of magnitude pQ . 

It should have been obvious from Fig. 1 that the decelerating force strength 
should be higher in absolute value than the accelerating one. 

However, on one hand, there is no known reason implying that a "social 
force" be instantaneous and have no duration. Moreover there is no reason that 



9 



the two growth rates, before and after an impulse, be identical, - the more so if 
the pulse has a finite duration 9. 

A different adaptation of the ideas in [1], on the evolution of competing 
entities, economic or sociologic ones, occurs if one writes 

dlr ^=k{l-X) + {a/e){l-X) (13) 



instead of Eq.(lO). In other words, one is (mathematically) letting Montroll's 
a to be A-dependent over the time interval [r; r + 9]. However, the emphasis 
differs much from [1]: rather than modifying the (Malhtus) X term, one adapts 
the (Verhulst) (1 — X) term to (economic or social) constraints. 
Mathematically, on the same footing as Eq.([i3|, one can write 



^ =kX(l-X) + (a/e)X(l-X). (14) 

Readily, the rate before the pulse is k but is k + a/9 after the "pulse" 
application. This "second rate" depends on the pulse strength and some time 
duration. 

Furthermore, the idea can be easily carried over within the Gompertz frame- 
work. It is "sufficient" to replace (1 — X) by ~ —ln(k/y); see Eq.([6]). A double 
log plot for N/Nm is shown in an appropriate replot of Fig. 2, i.e. Fig. 4, - the 
best fit equations being written in the figure. The fit is very precise. The lines 
are not parallel anymore. From these, one can deduce 9\ — 34 (months) and 
92 = 135 (months), i.e. ~ 3 and 11 years respectively. 

One can briefly elaborate on such observations. In fact, the cult present hi- 
erarchy interpretation and mine go along historical lines. The first acceleration, 
ca. 1914 can be historically connected to the first World War. 

After the war, income and housing taxes were implemented in Belgium, but 
social organizations were partially screened from such taxes. It took more than 
four years for a law on Organizations of Public Utility to be voted upon, - in 
1922. The parliamentary delays have in fine been much decelerating the temple 
construction and consecration process, such as introducing a remarkable time 
gap between the 16-th and 17-th temple consecration. After 1923, a smooth 
growth could resume, as seen in Figs. 1-2, with a lower but nevertheless appre- 
ciable rate. A remark should be made here: observe that there is not much 
effect ca. 1929, i.e. at the most famous financial crash time before the "Great 
Depression" . 

I wish to emphasize what seems to be a (financial) paradox in such a social 
force concept: the sect adepts were likely under much hardship due to the 1-st 
World War. Nevertheless, in need of some intra-community social support, the 
adepts donated quite an amount of money to build temples. Very interestingly, 
it has been found in a modern context that a similar situation, i.e. giving 
more at bad times, has just been occurring in Ireland [25]. This is a noticeable 
example of some sort of true altruism, in contrast to results from fund-raising 
campaigns 29 . 
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Figure 4: Gompertz plot : Log(Log) variation of the number of temples (N) 
raised in Belgium as a function of the number of months, starting from the raise 
of the first temple in Jemeppe-sur-Meuse, indicating the effect of a "social force" 
influencing a variation in growth rates, starting at t <~ 150. 
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4 Conclusions 



As an initial remark for a conclusion, let it be noted that both Verhulst and 
Gompertz first order differential equations for the evolution of a population are 
peculiar approximations of a more general one, i.e. 



d "=rN" 
dt 



y M J 



(15) 



where a, b, c can be requested to be positive real numbers; M is again the 
carrying capacity. Typical features on the relative and maximum growth rates, 
asymptotic values, inflection point characteristics are easily derived |16j . Intro- 



ducing x = (N/M) b , one can rewrite Eq.(15l as 
dx b r 



dt M 1 



x 1+ M\l-x] c . (16) 



Obviously, Verhulst equation is Eq.(15l or Eq.(16) when a=b=c=l. Gompertz 



growth equation is obtained for a=c=l, and b — >0. 

To introduce the Montroll's form, before applying a pulse force, one rewrites 
the above equation as 



d ln(x) 



dt M 1 



x 



(17) 



Thus, as long as the "social force" term has a structure similar to that of the 
right hand side, i.e. a power of x or of 1 — x, the integration can be easily 
performed., for example, in terms of incomplete Beta functions. This idea can 
be carried over within the Gompertz framework. It is "sufficient" to replace 
(1 - x) by ~ —ln(k/y)\ see Eq.(|6]). 

In summary, this study differs from that of Montroll who claimed that the 
evolution of some product results from a competition effect between two pro- 
cesses. He imposes a change in the x (or V)term, i.e. the already acquired 
content. In fact, Montroll's adaptation of Verhulst equation is more like an eco- 
nomic field competition description rather than a "social force". In the present 
considerations, the new term, called by analogy with Montroll's paper, a "social 
force" is attributed to a change in the limited need/capacity, in Verhulst or 
Gompertz equation, i.e. with emphasis on the (1 — X) term. The change in 
evolution is rather due to an external field, through the coefficient in front of 
(I-*)- 

In the studied example, the gap between the first and second regime has no 
doubt some extrinsic effect, i.e. an expected parliamentary law inducing a tax 
free system for the community, thereby hindering the temple constructions for a 
while. Thus, it can be pretended that one can be presenting two different views 
of a "social force": one resulting from a competitive interaction and another 
from an external field. 

Furthermore, note that the formal introduction of an extra (1 — X) term 
contribution rather than an X term based "pulse", as in Montroll's paper, not 
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only allows for greater flexibility in the data fit, but also in the analysis and in- 
terpretation of the resulting graphs. Indeed, the (1 — X) term adaptation allows 
for smoothness in the evolution, over some time interval, rather than requesting 
a jump in the evolution, with two constant and identical rates. Furthermore, it 
seems, "sociologically interesting", that the X term, like in Malthus approach 
bears upon what exists already, but the (1 — X) term, like in the Verhulst 
approach, bears upon what is missing. 

In conclusion, based on a relatively simple analysis of reliable data indicating 
adhesion to a growing sect, i.e. the construction of temples, a social force effect, 
in the physics sense introduced by Montroll, can be observed. However, a dif- 
ferent point of view can be taken. Moreover, comparing Verhulst and Gompertz 
law, it has been found that although both growth laws are based on different 
empirical considerations, the laws are found to be quite complementary From 
a sociological point of view of such complex systems, the " model" indicates that 
such communities are markedly influenced by external considerations ("external 
fields/forces"), beside intrinsic "religious" goals. 
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